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Waves and Jets in a Viscous Liquid, 

By A. B. Basset, M. A., F. R. S. 



1. In an able article which appeared in the ninth volume of this Journal, 
Prof. Greenhill discussed at considerable length the principal cases of wave- 
motion in a frictionless liquid which have hitherto been solved. In the present 
paper, I propose to consider certain problems of a similar character when the 
viscosity of the liquid is taken into account. 

The motion of .deep-sea waves is worked out in .Chapter XXIII of my 
Treatise on Hydrodynamics, from which it appears that when the viscosity is 
small, as is the case with water, the velocity of propagation is the same as that 
of a frictionless liquid, and that the only effect of viscosity is to introduce a 
modulus of delay whose value is equal to tf/Sn^v, where 3, is the wave-length 
and v the kinematic coefficient of viscosity. But when the viscosity is large 
the motion is of a totally different character, the time factor appearing in the 
form of an exponential quantity s u where Jc is negative. The result which is 
given in §521 of my book is wrong; for it can be shown, and will hereafter 
appear from the solution of another problem, that the correct value* of h is 
— gj 2rnv , whence the form of the free surface at time t after the commence- 
ment of the motion is 

£ = {a cos mx + b sin mx) e~2mv, 

where Injm is the wave-length, and a and b are constants depending upon the 
initial motion. The difference between the two kinds of motion may be illus- 
trated by allowing a jet of air to play for a short time upon the surface of water 
and upon treacle, and watching the effect produced. 

The motion of waves in a viscous liquid of finite depth is worked out in 
§523 of my Hydrodynamics. The solution is more complicated, but leads to 
similar results. 

* See also Tait, Proc. Eoy. Soc. Edinburgh, 1890, p. 110. 
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The motion of waves at the surface of separation of two frictionless liquids 
which are moving with independent velocities was investigated by G-reenhill* in 
1878, and his solution has been discussed by Lord Rayleighf with reference to 
the question of stability, but the corresponding solution for two viscous liquids 
would entail difficulties of a rather formidable character, owing to the fact that 
it would be necessary to retain all quadratic terms upon which the undisturbed 
motion depends. If, however, one of the liquids is frictionless whilst the other 
is viscous, the solution can be arrived at without difficulty. 

It is always a great advantage when the conditions of a mathematical 
problem can be expressed in terms of a single function, which satisfies a certain 
partial differential equation together with given boundary conditions. In most 
problems of interest relating to wave-motion in frictionless liquids this can be 
affected by means of the velocity potential ; but when the liquid is viscous, a 
velocity potential cannot as a rule exist, owing to the fact that the motion 
almost always involves molecular rotation. Of course any solution giving a 
possible irrotational motion must necessarily satisfy the equations of motion of 
a viscous liquid; but the solution will not satisfy the boundary conditions 
except in the single case in which the liquid moves like a rigid body having a 
motion of translation alone. Moreover, viscous motion involves a conversion of 
energy into heat. If, however, the motion of a viscous liquid is in two-dimen- 
sions, or is symmetrical with respect to an axis, the former can always be 
expressed by means of the current function ^ ; and this function will be used 
throughout the present paper. 

Waves at the Surface of Separation of Two Liquids. 

2. Let the origin be some point in the undisturbed surface of separation ; 
let the axis of x be measured in the direction of propagation of the waves, and 
the axis of z vertically upwards. Let the upper liquid be viscous and moving 
with constant velocity U lt whilst the lower is frictionless and moving with con- 
stant velocity TJ % . Let (i be the viscosity, and let the suffixes 1 and 2 refer to 
the upper and lower liquids respectively. 

Earnshaw's current function ^ satisfies the equations 

(^v , -w)v^ = (»^ + »'l)^ W 

•Cambridge Mathematical Tripos Examination. 
t Proo. London Math. Soc, Vols. X, XI and XIX. 
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where 

u = d4//dz, w = — d^/dx, (2) 

and v is the kinematic coefficient of viscosity which is equal to p/p. 

Since the liquid possesses an independent velocity TJ X , the solution of (1) 
must be of the form 

^ = 4/+Z7i 2 , 

where ty depends on the disturbed motion. Substituting in (1) and neglecting 
all quadratic terms which depend solely on the disturbed motion, we obtain 

(y-w-^y^o. w 

To obtain a solution suitable for representing wave-motion we assume that 
x and t enter in the form of the factor e imx + u j where 27t/m is the wave-length 
and k is a quantity to be determined. Substituting in (3) and putting 

rf — rrf + ih+imU^/v, (4) 

we obtain 

where $' now denotes a function of z alone. 
The solution of (5) is 



where % u % 2 satisfy the equations 

l -m 2 YA = 0,.} 

(6) 



&7h _•„ _ n 1 



3. We shall consider the case in which the depths of the two liquids are so 
large in comparison with the wave-length that the former may be treated as 
infinite ; under these circumstances the proper solutions of (6) will be 

Zl = As-™, %i = Or-, 

so that 

4 X = (Ae— + Ce— *) s imx + u + TJ x z . (7) 



(8) 
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The hydrodynamical equations for the pressure are 
du , _. du 1 dpi , „ 

dw . TT dw 1 dp, , „ 

-&•+*■& = — £±-v + «?"'> 

from these combined with (2), (3) and (6) we readily obtain 

Lm{h + imU l ) Zl = —-^- + 9, 

whence 

%>!= const. — gp lZ — ipiiJc + im^) Ae- m *+ im *+ u . (9) 

Since the lower liquid is supposed to be frictionless, the function t^g satisfies 
the equation 

whence 

^3 = .5e m *+ ' ma,+K + U 2 z (10) 

and 

jp* = const. — go % z + tp 8 (& + tm C7 2 ) Be mg+imx+M . (11) 

Having obtained the equations of motion and the pressure, we must now 
consider the boundary conditions which hold good at the surface of separation. 
Let the equation of this surface be 

£ — ae imx+ * t =F(Z,x,t) = 0, (12) 

then the kinematical condition which holds good at every bounding surface is 

dt ■" dz dx dx d£ 
Applying this to the upper liquid we get 

a= _ ^(A + c) (14) 

k-\- imUi 



Applying (13) to the lower liquid we get 

LtnB 



a = 



(15) 



The dynamical conditions which hold good at the boundary are continuity 
of stress ; and since a frictionless liquid is incapable of supporting any tangential 
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stress or shear, it follows that the tangential stress parallel to the plane 2=0 
must vanish at the surface of separation. This condition gives 

du t dw . 

dz dx 

£-£=<>. <"> 

when 2=0. 

To obtain the second condition, consider a film of the viscous liquid 
bounded on the lower side by the surface of separation and on the upper side 
by a surface indefinitely near the former. The lower surface of the film will be 
subjected to the pressure j> 2 of the frictionless liquid, and the upper to a traction 
B due to the action of the viscous liquid in contact with it ; hence if T be the 
surface tension between the two liquids, we get 

p , + B = j=-T§, (17) 

since the curvature p -1 of the film is supposed to be very small. 
Substituting from (7) in (16) we obtain 

2Atn? + (a 2 + O = . (18) 

The value of B is 

B =-.ft + ^-ar = -!*-%*' ass 

= 9?x a + 'Pi (& + tm Ui) A + 2ip 1 mv [Am •+■ Oa) 

by (9); whence by (11) and (12), (17) becomes 

9 (pi — P2) a + tpi (k + t m Ui + 2mV) A 

+ 2iprnavC +(,&(& + irnUz) B— 7 T wi 2 a = 0. (19) 

The elimination of a, A, and B between (14), (15), (18) and (19) leads 
to an equation for determining k. 

From (4), (14) and (18) we obtain 

t A = — a(h-\- itnUi + 2<m % v)\m, iG= 2amv; 
whence the result of the elimination is 
m 9{?\ — 9%) — pi(k + imU 1 + 2m % vY + Apipfav* — p 2 (& + im,V 2 f — Trn 3 = 0. (20) 
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Equation (20) combined with (4) gives a biquadratic equation for deter- 
mining k, and we shall first perform the elimination. 

Let x = h + cmUi (21) 

and 

1= (pi + p 2 ) x % +2x{ 1o x m\ + LO % m{U % — 0i) } 

+ mg( Pz -p 1 ) + Tm s -o s m\U z -U l f, (22) 
then (20) becomes 

J= 4p 1 mV (a — in)) 
also (4) may be written 

x = v (a 2 — rn % ) , 
whence 

P + 8 JpxmV — 1 6pjm Vsc = . (23) 

This is the biquadratic equation which we have to discuss. 

Oil Floating on Water. 

4. The viscosity of water is equal to the tangential stress of .014 dynes 
per square centimetre ; and since a dyne is approximately equal to the attrac- 
tion of the earth at its surface upon a cubic millimetre of water,* the viscosity 
of water is a very small quantity. On the other hand the viscosity of olive oil 
is about 3.25 dynes per square centimetre, and is therefore about 232 times as 
great as that of water. If, therefore, a thick layer of oil is floating on water, 
the former may be treated as a highly viscous liquid and the latter as a friction- 
less one, and consequently the proper solution for our purpose will be one in 
which the viscosity is supposed to be a large quantity. 

Assuming for trial that x varies as v~\ it follows that the first term and 
also the imaginary term in i" are negligible ; accordingly (23) reduces to 

J= 2p 1 wi ! W, 
which gives, since vpi = [t, 

(i). Let Z7i = U % = , so that the liquids are initially at rest ; then, since 
P2>pi) & i s a negative quantity and wave-motion in the proper sense of the 
word does not exist. The amplitude of the initial disturbance gradually dimin- 
ishes with the time and the disturbed motion ultimately dies away. 

* More accurately the 981st part of a cubic centimetre. 
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When the wave-length is large, m or 27t/Jl is small, in which case the term 
depending on gravity is the most important ; on the other hand when the wave- 
length is small, the term depending on the surface tension is the most important. 
It therefore follows that the damping effect of gravity is most important in the 
case of long waves, and that of surface tension in the case of short ones. 

We have supposed that the upper liquid is oil and the lower one water, so 
that p 2 >• pi ; if, however, ^ > p 3 so that the upper liquid is the heaviest, the 
value of Jc will be positive unless 

IW^gfa — p,), 

and consequently unless this inequality is satisfied the equilibrium will be 
unstable. If, however, the difference of densities is very small and the waves 
are sufficiently short, it will be possible for this inequality to be satisfied, and 
we shall have the curious case of a heavier liquid floating on a lighter one in 
stable equilibrium. The stability will, however, be confined to disturbances of 
very short wave-length. 

When the wave-length is such that the numerator of (24) vanishes, a critical 
case arises. This case is best discussed by reverting to the original differential 
equation for k. It, will be found that one root is zero and that another is real 
and negative, so that there must be one stable mode of motion, but this case is 
not of any particular interest. 

(ii). When p 2 ^p! and TJ X , U % are not zero, it follows that the real part of h 
will be positive provided the term ( Z7 2 — U-tf is large enough. From this it 
follows that a high relative velocity, which includes the case of motion in oppo- 
site directions, tends to produce instability. This has been verified experi- 
mentally by Prof. Osborne Reynolds.* The exact nature of the instability will 
depend upon whether the two liquids are capable of permanently mixing. If 
they are, the subsequent motion will at first be of a highly turbulent character, 
which will continue until the liquids have become thoroughly mixed together ; 
but if they are incapable of remaining permanently mixed, a layer of turbulent 
liquid will exist in the neighborhood of the surface of separation, resembling a 
mixture of oil and water which has been thoroughly well shaken together. If 
the liquids be left to themselves, the damping effect of viscosity will gradually 
cause the turbulent motion to subside and the ultimate state will be one in which 
all relative motion between the parts of the viscous liquid will have disappeared. 

*Phil. Trans., 1883. 
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It must, however, be recollected that as the lower liquid is supposed to be fric- 
tionless, it is not necessary that the tangential velocities on either side of the 
surface of separation should ultimately become equal. 

The critical case in which the numerator of (24) vanishes would require 
separate treatment. 

5. We shall now consider the case in which the viscosity is very small, and 
we shall suppose that the lower liquid is viscous. This renders it necessary to 
reverse the sign of g. Since v is very small, the last two terms of (23) may be 
neglected, and the equation for k is 

(pi + p.) & + 2* { fym'v + tp 2 m ( U % — J7i) f 

+ *ng(pi-p>) + Tm s -o 2 ?n*(U 2 - Utf = 0. (25) 

If the sum of the last three terms is positive, it follows that the roots of 
this equation must be of the form 

a -f t/3, y{a — i/?), 

where y is a positive quantity ; also 

4p 1 mV = — a (I + y) , 

which shows that a must be negative ; hence the real part of h will be negative, 
and the motion will be stable, provided 

9 (pi - p») + Tm* > o % m ( Z7 2 - U,f. 

If this inequality is not satisfied y will be negative, and therefore one root of 
h must be positive, so that the motion will be unstable. 

When both liquids are at rest and p x > p 3 , the real part of h is equal to 
— 2p 1 wr ! j;/(p 1 + pa), which shows that the amplitude diminishes with the time, 
and the wave-motion gradually dies away. 

Motion of a Sheet of Liquid. 

6. In §418 of my book, a short investigation is given respecting the motion 
of a thin sheet of frictionless liquid. We shall now discuss the same problem 
when the liquid is viscous, and surface tension is taken into account; but in 
order to avoid unnecessarily complicating the analysis, the effect of gravity will 
be neglected. The conditions of the problem may be approximately realized in 
practice by supposing liquid to escape from a long and narrow horizontal slit in 
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the side of a cistern, and confining our attention to that portion of the sheet 
which is nearly horizontal. 

Let U be the undisturbed velocity of the sheet, 2l its thickness. We shall 
also suppose that the disturbed motion is such, that the sinuosities are sym- 
metrical with respect to the middle surface of the sheet. 

The value of -fy will be 

where % lt % t are determined by (6) of §2. Since we suppose that the motion is 

symmetrical with respect to the plane 3 = 0, d-^/dz must be an even function 

of a, and therefore 

$ = (A sinh me + sinh ass) g"»*+" + Uz, (1) 

whence from (8) of §3 we obtain 

p = ip(k + imU)A cosh nize im * +kt (2) 

Let the equation of the free surface of the upper side of the sheet be 

f— I — as"- m!e+kt — F— 0, 

then applying (12) of §3, we get 

mi {A sinh»i?+ (7sinhaZ) , oS 

a = J j—. ff -. (3) 

The condition of zero tangential stress gives 

lArr? sinh ml + G{a? + m z ) sinh a? = . (4) 

Again, omitting the exponential factor, 

R=-p-2 v?mz 

= ip(h + tmU)A cosh ml — 2ipmv {Am cosh ml + Ca cosh at) . 
The surface tension condition is 



B = -~= T 



a ' 



dx' 
which becomes 

io (Jc + i m U+ 2m V) A cosh ml + 2ipvmaC cosh al — Tm z a. (5) 

On applying the same conditions to the lower side of the jet, where 
z = — I — ae imx+M , it will be found that they lead to exactly the same equations ; 

hence on eliminating a , A and G we obtain 

Tm 3 

(7e + LmU+ 2m 2 vf coth ml — 4m 3 v 3 a coth al -\ = . (6) 

r 
This equation combined with (4) of §2 determines the values of h. 

14 
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When v is very small, a 2 = (h + tmU)/v approximately; a is therefore a 
large quantity, and coth aZ= 1. Also av % is very small, so that the second term 
of (6) may be neglected. Under these circumstances we obtain 

iTm 3 
h= — 2m?v — imU±:iJ t&nh ml. (7) 

From this equation we see that the real part of h, which is proportional to the 
viscosity, is negative, and therefore the motion is stable, from which it follows 
that the tendency of viscosity is in the direction of stability. 

We also see that the effect of surface tension is to introduce an imaginary 
term into the value of 7c, which consequently does not affect the stability. In 
the case of a circular jet, it is known that surface tension produces instability 
when the disturbance consists of short waves. It is remarkable that surface ten- 
sion should produce such dissimilar effects ; and it is important to notice that 
the effects, produced by certain causes upon a liquid whose motion is in two 
dimensions, are often totally different from those produced by the same causes 
when the motion is in three dimensions. The former kind of motion is generally 
the simplest from a mathematical point of view, whilst the latter kind more 
often occurs in practical and experimental applications. Great caution must 
therefore be observed in drawing inferences from theoretical results derived 
from the solutions of problems in two-dimensional motion and in applying them 
to the explanation of experimental phenomena in which the motion is in three 
dimensions. 

When v is absolutely zero, equation (7) gives the value of k for a "friction- 
less liquid. This result can be verified by an independent investigation. 

Equation (6) may be written in the form 

(J+iw Uf coth ml + 4m?v (h -f imU) coth ml 

Tm 3 
-f 4m 3 v i (m coth ml — a coth al) -\ = 0. 

When v is very large, a is nearly equal to m, so that we may write 

a = m + (& + tmU)/2mv 

approximately ; whence, by Taylor's theorem, 

la coth al = Im coth ml + - — ~ — (coth ml — ml cosech 2 ml) , 

and the equation for h becomes 

(h + i m TJ) % coth ml + 2m 3 lv (k+ imU) cosech 2 ml + Tm 3 /p = , 
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which gives 7 , _ SPsinh 2 ml 

k + cmU= ^— 

approximately. From this result we see that when the viscosity is large, its 
effect, combined with that of the surface tension, is to cause the disturbed motion 
to die away. 

Motion of a Cylindrical Jet. 

7. The motion of a cylindrical jet of frictionless liquid was discussed by 
Lord Rayleigh about ten years ago f and in two recent papersf he has returned 
to the subject and considered the case of a highly viscous liquid. I now propose 
to. give an investigation of a far more general character, in which the velocity, 
viscosity, surface tension and also the influence of the surrounding air are taken 
into account. I shall in addition suppose the jet to be electrified. 

When a fine jet is allowed to escape from a nozzle so as to describe an 
approximately parabolic path, it is observed that at a certain point the jet 
becomes disintegrated into drops. If, however, a metal ball, which is so feebly 
electrified that it scarcely causes the leaves of a gold-leaf electrometer to diverge 
to any perceptible extent, is brought into the neighborhood of the jet, the disin- 
tegration ceases. It is, moreover, necessary for the electrification to be feeble, 
for a powerful charge causes the disintegration to proceed more rapidly. I 
believe that Lord Rayleigh considers that this phenomenon cannot be completely 
explained by ordinary mechanical methods, but is mainly due to the fact that 
the electrification causes the drops to coalesce instead of rebounding upon col- 
lision. At the same time very little complication is introduced into the work 
by supposing the jet to be electrified, and the results are instructive, as they give 
information as to the effect which a charge is likely to produce. We shall find 
that the charge tends to produce stability or instability according as the ratio of 
the circumference of the jet to the wave-length of the disturbance is less or 
greater than .6 . 

8. We shall suppose throughout that the jet consists of a straight cylindrical 
column of liquid, and that the disturbance is symmetrical with respect to the 
axis. Taking the axis of the jet as the direction of the axis of z, and measuring 

* Proc. Lond. Math. Soc, Vols. X, XI and XIX. t Phil. Mag., Aug., 1892. 
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r perpendicularly to it, the motion can be determined by Stokes' current function, 
which, in the case of a viscous liquid, satisfies the equation 

where 

D= — + — — — — 

and 1 ^ 1 <ft// /0 \ 

r dz ' r dr 

When the motion is undisturbed, w = TP, where IP is the constant velocity 
of the jet j whence, from the second of (2), 

which is readily seen to be a solution of (1) . 

Let V be the portion of i|/ which depends upon the disturbed motion, so that 

Substituting in (1) and neglecting all quadratic terms which depend solely upon 
the disturbed motion, we obtain 

In order to obtain a solution which represents wave-motion, we assume 
that z and t enter in the form of the factor e im * +u , whence if 

a 2 =m 2 + (k + imW)/v (4) 

(3) becomes 

( ' $ 1 d aV d 2 Id ,\ ,, n . . 

Kd^-V-W- m )W-ir^r-o?)^ = o, (5) 

the solution of which is 

•V = Zi + Z»> 
where 

^&__L^?_ aV -o f 

In the first of (6) write % t = rty , and the equation for <p becomes 



(6) 
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The solution of this equation is 

$ = AIi (mr) -f BK X (mr) , 

where I lt K x are the two associated Bessel's functions of order unity, and are 
connected with the functions of order zero by the equations 

/„' = />, Kl = K u (7) 

where the accent denotes differentiation with respect to the argument mr. These 
functions are discussed in §§266-269 of my Hydrodynamics, and also in Chapter 
XIX of Lord Rayleigh's treatise on Sound. 

The function I (x) is expressible in the form of the series 

Ct* Or Ou 

-*o( a 7 = 1 H~ IjE" + 22 42 "i~ 2 2 4 2 6 a "^ " " " ' ' 

whence r . s x , x 3 , x B , /f o. 

Ii(x) = -3- + ^4 + 2X4175 + (9) 

whilst the function ^" (a;) can be expressed by means of either of the series 

K (x) = (y + \og^x)I (x) ^-S 1 — -^-^S 2 —^^-^S i — , (10) 

E *( x ) = -(te) a J 1 " T78x+T^(8xT- 1.2.3.(8*)' ~\ (11) 

where y is Euler's constant, which is equal to .5772 .... , and 

Since the / functions are finite when x = , they are suitable for space 
inside a cylinder ; on the other hand the K functions are infinite when x = 
and zero when a5 = », and are therefore suitable for space outside. It therefore 
follows that the proper value of <i> is 

+ = r{AI 1 (mr) +CI 1 (ar)\ e im ' +7 * +$Wr*. (12) 

In taking account of the influence of the surrounding atmosphere we shall 
neglect its compressibility and viscosity ; whence, if ^i be the current function, 

D^ = 

and 

^ 1 = BrK 1 (mr) «""+*. (13) 
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The hydrodynamical equations for the pressure are 

du du 1 dp , / , u \ 

-dt +w ^ = -yd7 +v (y u --wj' 

dw dw 1 dp „ 

The first of these equations may be written 

tm 1 dp t / d? 1 d 1 \ i>' 



~~ o dr ' r Vrfr 2 r <2r */' 



P 

whence 1 dp irnv , „ „. 

p a/r r v ' ~ 

and _£>= ip(& + tmW)-zLr o (jnr) + const., (14) 

which can also be shown to satisfy the second equation. 

In the same way the pressure of the air can be shown to be 

p x = ipihBK (mr) -f- const. (15) 

9. We must now consider the boundary conditions. 
Let the equation of the surface of separation be 

a + Gs' mi + kt — r = F=0, (ISA) 

then the kinematical condition at this surface is 

i^__L W §E. i J_ _# dH — o (16) 

dt r dz dr r dr dz 

Applying this to the jet we obtain 

— tffl Mji ( ma ) "t~ ^ ( aa ) I ( 17 i 

C ~~ k + imW ' l ; 



Applying (16) to the surrounding atmosphere we get 

imBKi (ma) 
e - _ . 



(18) 
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The dynamical 


condition of 

du 
dz 


zero tangenti 

. dw 
+ 37 = °' 


al stress 


gives 


or 




d?V 
dr 1 


l_d4> 
r dr 


«jV' = 


0, 


or 




2m? Zl 


+ (a 2 + m 2 ) 


£2 = 0, 




and therefore 











ton 9 Al,. M + (a 3 + m 8 ) GI X (aa) = . (19) 

To find the remaining boundary condition, let V be the potential of the 
induced charge; then the resultant force which acts outward on. an element of 
the jet is 

__1_ dV_ 1 /dV\* 
2 dr ~ 8?t \ dr J ' 

also the portion due to surface tension measured outwards is 

\ p x p 2 a J 

where p x is the radius of curvature of the meridian section of the disturbed sur- 
face, and o % the corresponding quantity for the perpendicular section. If there- 
fore we consider the equilibrium of a film of the disturbed surface, the required 
condition is easily seen to be 

NOW D — ^JL.O du 



d /^'\ 
= -p-2cpmv w (jr) 



Substituting the value of p from (14) and the value of 4' fr<> m (12). an d 
recollecting that 



I[(x) = I(/=-^- + I , 



we obtain 



P=z — ip(k+ tmW + 2m 2 v) AI (ma) + 2ipmva~ 1 AI 1 (ma) 

— 2ipmva~ l \aal (aa) — I x (aa)} + const., (21) 

the exponential term being omitted for brevity. 
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Again, from (15 A), 

JL _ _ ®L — 2 

p x dz % ~ ' 

1_ l_ _ 1 _1_ _ _c_ 

p 3 a a + c a « 2 ' 

whence the surface tension term becomes 

f(mV-l)c/a ! . (22) 



The value of F is 



7= 



ffg (wir) ^ 



2^ log r + " v s 
& ^ K (ma) 



£' 



where i? is the charge on the undisturbed jet per unit of length, and F is a con- 
stant to be determined from the condition that V is constant when r is given 
by (15 A). Dropping the exponential factor, this condition gives 

F= 2Ec/a, 
whence at the surface 

_dV _ 2E 2Ec maK[(ma) 
dr ~ a + c a 2 K (ma) ' 
and therefore 

/dV\*__4E* ±E*c< maK^ma) ) 

KdFj-T? a ir \ l+ K (ma) ]' W 

Substituting from (21), (15), (22) and (23) in (20), and equating to zero the 
coefficients of the omitted exponential term, we finally obtain 

— ip(k + imW + 2mV) AI (ma) + 2ipmva~ 1 AI 1 (ma) 

— 2ipmva~ l O \ aal (aa) — Ii (aa) \ + ipJcBK (ma) 

The elimination of the four constants A, O, B and c between (17), (18), 
(19) and (24) furnishes the required equation for determining h. The result of 
the elimination gives the following equation for k: 

(k + l mW+2m*vr^^ ) -^(Jc + LmW+2m*v) 



I (aa) 4wV _ ffptifo (ma) T(m z a % —1) 
^(aa) "■" a pmK x (ma) "*" pa 2 



+ 



E % f maK-t (ma) ) 



4npa' 
where p 3 denotes the density of the atmosphere. 



< maKAma) ) 

\ ^ K (ma) J v ' 
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10. We shall now discuss this equation when the viscosity is so small that 
powers of v above the first may be neglected. 

Let q = h + imW, x=ina; then from (4) it follows that a = (q/v) i 
approximately, provided the wave-length is not so small that in 2 is comparable 
with q/v.- Equation (25) may now be written 

q% (x) 2v_ j 2xl (x) "I 4mV j x J (x) aal (aa) \ 

xl,{x) + a? \ I x {x) 1 \ q+ a* } I^x) 7,H J 

_^w + n<f 7 D + j f 1+ -gwi = , (26) 

oxK x (x) pa 3 ' 47tpa s [ K (x) \ v ; 

Now if a be very large, 

from which it follows that the limit of Lj/I-l when a = o° is unity. The third 
term of (26) accordingly does not involve any lower power of v than v i , and may 
therefore be neglected. 

11. We shall now consider the effect of the last three terms of (26) sepa- 
rately. 

(i). When Zand E are zero, (26) may be written in the form 

& 8 f p/o(g) Pi-g"o(*) ) , 2 y ftf 2sZ (a!) ) 2 t W t (x) W%(x) _ . 

p (^(a;) - xZi(a;)J "+" a 2 ( ^(a?) j "*" xl^x) xl^x) ~ V ' K '' 

Now 7 and 7 X are always positive, but the quantity fi^/iTi is always nega- 
tive, whence the coefficient of Ti must be positive ; also from the series for I , I± 
it follows that the second term of (27) is also positive, whilst the last term is 
negative ; whence the roots of (27) must be of the form p (1 + t/?) , y (1 — i(3) , 
where the two quantities p and y are of opposite sign. Hence the real part of 
one of the values of h is positive, and therefore the motion is unstable. 

(ii). The fifth term of (26) depends upon the surface tension. If x is greater 
than unity, this term is positive, and consequently the real part of Tc must be 
negative and the motion is stable. If the liquid is supposed to be frictionless, 
so that v is absolutely zero, we fall back on the case discussed by Lord Rayleigh, 
in which it is shown that the motion of a jet will be stable or unstable according 
as the ratio of its circumference to the wave-length of the disturbance is greater 
or less than unity. 

When x<C 1, the fifth term of (26) is negative, and consequently one value 
15 
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of h must be positive ; but in order to examine the case of values of x which are 
very slightly less than unity, it would be necessary to take into account the 
third term of (26) which has been neglected. This would be troublesome, but 
there is little doubt that its effect would be found to be that motion for which x 
is slightly less than unity would be stable ; in other words, the effect of viscosity 
would be to render the liquid less unstable than if it were frictionless. 

(iii). We must next consider the electrical term. 

The functions K Q (x) and xK x (x) have been tabulated by Lord Rayleigh* 
From this table it appears that K (x) is equal to — oo, — .9244, and — .7774 
for x = oo, .5 and .6 respectively, whilst the corresponding values of xK x (x) are 
1 , .8283 and .7817. From these figures we see that xK x jK§ is zero when x = 0, 
and becomes equal to —1 for a value of x which is slightly less than .6. Ac- 
cordingly for disturbances whose wave-length is large in comparison with the 
circumference of the jet, the charge tends to produce stability, but for disturb- 
ances of short wave-length the charge tends to produce instability. 

(iv). The effect of a very small viscosity considered by itself is to produce 
stability, for in this case one value of q is zero, whilst the other is 



2^ a ( - Ii(x) \ 
a* Y xl (x)y 



The conclusions that we draw from the preceding analysis are the following. 
Viscosity tends to produce stability ; the influence of the surrounding air tends 
to produce instability ; the surface tension tends to produce stability or insta- 
bility, according as the ratio of the circumference of the jet to the wave-length 
is greater or less than unity; the electrical charge tends to produce stability or 
instability according as the above ratio is less or greater than .6. 

These results do not entirely explain why it is that a very slight charge 
produces stability, and it must be admitted that the phenomena which are 
actually observed depend upon some cause of which account has not been taken, 
but the results do appear to give a satisfactory explanation of why it is that a 
large charge produces instability, for the charge produces instability for wave- 
lengths for which x >• .6, and when it is sufficiently great its effect will over- 
come the steadying effects of viscosity and also of surface tension. 

The effect of a large viscosity has been discussed by Lord Rayleigh, so that 
it is unnecessary to say anything upon this point. 

•Phil. Mag., Aug., 1892, p. 179. His <j>(z) and f (x) correspond to my K (x) and K^x) . 



